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Mesonic quasinormal modes of the Sakai-Sugimoto model at high temperature 

Nick Evans 1 - 2 and Ed Threlfall 1 

1 School of Physics and Astronomy, University of Southampton, Southampton, SO 17 IB J, UK 

2 Newton Institute for Mathematical Sciences, 20 Clarkson Road, Cambridge, CB3 OEH, UK 

evans@phys.soton.ac.uk, ejt@phys.soton.ac.uk 

We examine the mesonic thermal spectrum of the Sakai-Sugimoto model of holographic QCD by 
finding the quasinormal frequencies of the supergravity dual. If flavour is added using D8-D8 branes 
there exist embeddings where the D-brane worldvolume contains a black hole. For these embeddings 
(the high-temperature phase of the Sakai-Sugimoto model) we determine the quasinormal spectra 
of scalar and vector mesons arising from the worldvolume DBI action of the D-brane. We stress 
the importance of a coordinate change that makes the in-falling quasinormal modes regular at the 
horizon allowing a simple numerical shooting technique. Finally we examine the effect of finite 
spatial momentum on quasinormal spectra. 



INTRODUCTION 

Gravity dual descriptions [3, 0, 0| of strongly coupled 
gauge theories with quarks have recently shed light on 
the physics of mesons and chiral symmetry breaking JiJ. 
There has also been considerable interest in studying the 
finite temperature behaviour of these systems. 

The simplest such dual is the near horizon geometry of 
a D3-D7 brane system which describes an M = 2 gauge 
theory. The D7 branes can be treated as probes || in the 
limit where the number of flavours is much less than the 
number of colours, -WV<C N c , or the full back reacted ge- 
ometry can be found @,[Bi[3| ■ The meson spectrum in the 
probe limit has been computed in Q . Finite temperature 
manifests itself as the presence of a black hole in the dual 
space-time [2j . In the infinite volume limit the black hole 
geometry is energetically preferred for any temperature 
greater than zero. The transition from an AdS space cor- 
responds to the analogue of the deconfinement transition 
in the pure glue gauge theory (which in a conformal the- 
ory occurs as soon as the dimensionful parameter, T, is 
introduced). A further first order phase transition has 
also been found in this system fiol . HH [tH . 13, 14 1 when 
the temperature passes through the scale of the mass of 
the mesonic bound states - this corresponds to when the 
horizon of the black hole grows to swallow the D7 probe 
in the interior of the space. This transition has an as- 
sociated small jump in the chiral condensate's value but 
the main physics of the transition appears to be the me- 
son fields melting into the thermal background. Once 
the D7 brane enters the horizon there are no longer nor- 
malizable fuctuations of the D7 brane that generate a 
discrete set of meson bound states. Instead there are 
quasi-normal modes of the black hole corresponding to 
fluctuations of the D7 brane that are pure in-falling at 
the horizon. These fluctuations correspond in the dual 
gauge theory to excitations of the plasma with a complex 
mass parameter - the excitations have both a mass and a 



decay time. The spectrum of these quasi-normal modes 
has been explicitly computed in poj . 

Another interesting model is the Sakai-Sugimoto model 
lil ljj which is based on a (wrapped) D4 D8 D8 system. 
It is a gravity dual of a non-supersymmetric gauge theory 
(which is four dimensional in the IR but five dimensional 
in the UV) that dynamically breaks a non-abelian chi- 
ral symmetry of its quark fields. The high temperature 
phase again corresponds to a transition to a black hole 
geometry. The transition occurs when the black hole's ra- 
dius becomes of order the wrapped circumference of the 
D4 brane which is also the parameter that determines 
the mass gap of the theory. This behaviour is more akin 
to what one would expect in QCD than that of the con- 
formal theory discussed above. 

Massless chiral quarks can be introduced by placing the 
probe D8 and D8 branes at anti-podal points on the cir- 
cle the D4 brane is wrapped on. In the near horizon 
limit of the D4 branes these D8 branes choose to join 
at the scale of the mass gap breaking the chiral symme- 
tries on their world volumes to the diagonal sub-group 
and generating a mass gap for the mesonic fluctuations 
of the D8s. When the geometry makes the transition at 
finite temperature to the black hole background the D8 
and D8 disconnect and instead lie straight and fall into 
the horizon [l9| . Chiral symmetry breaking is therefore 
restored along with deconfinement. 

There is a larger class of embeddings in which the D8 
and D8 join at a larger radius in the space so there is a 
bigger mass gap for the quarks. In [171 ] we have argued 
that these embeddings describe a quark mass in the the- 
ory although it has been also argued in the literature 
[ 1 81 ] that the chiral symmetry breaking scale is being en- 
hanced in these cases by higher dimension operators. The 
distinction is not important for what we discuss here - in 
these cases there is a further first order transition as the 
temperature (horizon) grows through the mass scale of 
the mesons. The transition is very much like that of the 



D3-D7 system in that the mesonic fluctuations of the D8 
branes are replaced by quasinormal modes of the black 
hole. The mesons of the theory have melted into the 
plasma. 

The mesonic fluctuations of the D8 branes above the 
phase transition have been studied in 2l|. Here we will 



concentrate on the very high temperature phase where 
the D8 branes lie straight and fall into the black hole hori- 
zon. We will explicitly compute the quasinormal mode 
spectrum corresponding to the scalar and vector mesons 
of the theory. 

As a prelude to this we compute the quasinormal spec- 
trum of a Klein-Gordon scalar living on the D8 brane 
worldvolume. We apply the idea of regularizing the co- 
ordinates for ingoing modes, which has previously been 
used in asymptotically-flat spacetimes (see for example 
[H) and in the context of AdS-CFT in [H. The result 
of this is that the ingoing mode is described by a regular 
Taylor series at the black hole horizon. We use this as 
the initial condition and obtain the quasinormal spectra 
by shooting out from the horizon. We wish to stress that 
this is a much cleaner numerical process than trying to 
match on to oscillating solutions at the horizon. We use 
the same method to examine the spectra of modes arising 
from the DBI action of an embedded D8 brane. We treat 
a scalar fluctuation of the brane in the geometry and a 
Lorentz vector arising from the Maxwell field on the D- 
brane. Finally we briefly discuss the effect of nonzero 
momentum on the spectra and extract the diffusion coef- 
ficient from the lowest quasinormal modes of longitudinal 
vector excitations in the small k 'hydrodynamic' limit. 



THE GEOMETRY 



The metric of the high temperature Sakai-Sugimoto 
model is 



ds 2 = (% 



dx\ 



dT 2 



Here f(u) 



-f(u)dt 2 

, (R\% ( du 2 
^Ku) \f(u) 

1 — (^j 1 ) 3 and the dilaton is e" 



u 2 dVL\ 



(1) 



The parameter ut, representing the position of the hori- 
zon in the geometry, gives the temperature in the dual 

i 

field theory by the relation T — This is the Hawk- 

R^ 

ing temperature of the black hole. 

Let us work in the dimensionless radial coordinate x = 
— and measure the Minkowski and r dimensions in units 



°^ v «r ' This corresponds to measuring frequencies and 
momenta in units cx T. 

The metric becomes 



R? Jut 



X 2 



+X 2 



f[x)dt 2 
( A 



dx 3 



dT 2 



(2) 



/O) 



x 2 dfll 



We will consider the spectrum of modes associated to a 
D8 brane. The D8 branes fill the space except for the 
r direction in which they live at a single value of r = 
To - this is the energetically preferred high temperature 
configuration [ijj]. The induced metric on the D8 brane 
worldvolume is just (T2|) with dr = 0. 



REGULAR HORIZON COORDINATES FOR 
QUASINORMAL MODES 



As a warm up we shall first consider linear fluctuations 
of a Klein-Gordon scalar restricted to the worldvolume 
of the D8 brane. This is not a physical mesonic state of 
the field theory but it exhibits similar quasinormal modes 
in the supergravity dual and we use it to illustrate our 
calculational technique. The 9D action describing the 
fluctuation is 



Ssf — — 



(3) 



Here the geometry is given by the induced metric on the 
D8 brane. 



The equation of motion for the fluctuation is 
1 



-gg^M = 



(4) 



Writing our the equation for a scalar fluctuation with spa- 
tial momentum k with respect to the plasma rest frame 
and zero 64 spin as $ oc e _ **' t+tfe ' a:3 one obtains the equa- 
tion 



■ X 4 



$ = 



(5) 



Here the prime indicates an x derivative. The large- a; 
asymptotic of the equation is =0 with solution 

$ ~ ci + C2 x~~^ . For a normalizable solution we clearly 
want the decaying power. 

Taking the near horizon limit (x — ► 1) one finds $ ~ {x — 
l) ±l ^ . Since the whole solution is exp(— iuit±if ln(x— 1)) 
the solutions are ingoing (— ) and outcoming (+) waves 



(I 



Klein-Gordon scalar 



lm((i>) -3 ; 

-4 ; 

-5 1 



II 




1 


± 2.9348 - 1.6032 i 


2 


± 4.3133 - 2.3760 i 


3 


± 5.6509 - 3.1367 i 


4 


± 6.9725 - 3.8928 i 


5 


± 8.3151 - 4.6531 i 
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FIG. 1: The lowest five k — Klein-Gordon scalar quasinormal frequencies in the complex u plane. 



which oscillate infinitely many times before reaching the 
horizon. 

In these coordinates it will be hard to find the full solu- 
tion numerically - one needs to shoot out from, or on to, 
a highly oscillatory solution near the horizon. A better 
way to proceed is to change coordinates to make the in- 
falling solution regular at the horizon so numerical meth- 
ods can be more easily used. In particular we will shift 
coordinates so 

t = h — a(x) (6) 
To make the infalling solution regular we require 



a(x) — — \n{x — 1) + 
3 



(7) 



where the additional terms are regular at the horizon. 
One way of satisfying this is to define 



1 



da 



dx x 3 — 1 



(8) 



We can now recompute the scalar equation of motion and 
we find 



x * /$ 



7 3 

2x±& + -x 4 $ 



(a/ -r)^$ = 

(11) 

We note that this is an equation with five singular points 
in the complex x-plane. The general solution of such 
equations is not given in terms of well-known functions, 
nor can one immediately apply the continued fraction 
method as used in [24| . Accordingly we use a purely 
numerical approach. 

There are two criteria for a good solution. Firstly the so- 
lution should be a purely ingoing wave at the black hole 
horizon since classically a black hole can absorb but not 
emit particles. Secondly the solution must be normaliz- 
able when integrated over the D-brane worldvolumc. 

The large- a; asymptotic of the equation is the same as 
we found in the usual 'Schwarzschild' coordinates and so 
we choose the decaying power which is normalizable as 



We then have 



dh = dt 



1 



dx 



x 3 f(x) 

leaving the metric in the new coordinates as 



At the horizon we will seek a solution in the form of a 



Frobenius series $ 



1. Now 



(9) substituting this into the differential equation yields the 
indicial equation 



ds 2 



-I! f(x)dh 2 + 2x zdhdx + x ^dx 2 



(da 



dT 2 



y/xdfl 



(10) 



The induced metric is the above with dr = 0. Note that 
the near horizon solution of ([5]) only depended on the 
powers of / in the function not the powers of x (which 
becomes one in the near horizon limit). These powers 
of / will turn out to be the same for all of the modes 
we consider below and so this change of coordinates will 
suffice to make all infalling modes we look at regular. 



2 



(12) 



The solution with s = is regular at the horizon and 
corresponds to a purely infalling solution. 

Using the regular Taylor series for the infalling solution 
as the initial condition we shoot out from the horizon. 
By requiring our solution to vanish as x — > oo we can 
find the quasinormal frequencies. They are displayed in 
Fig.l. 
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n 
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± 3.3896 - 1.9498 i 


2 


± 4.7560 - 2.7320 i 


3 


± 6.0900 - 3.4970 i 


4 


± 7.4110 - 4.2297 i 
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± 8.7234 - 5.0105 i 
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FIG. 2: The lowest five k — sigma meson quasinormal frequencies in the complex u plane. 



QUASINORMAL MODES FROM FLAVOUR 
BRANES 

We will now consider the physical quasinormal modes of 
D8 branes, coming from the DBI action. 

The presence of the brane corresponds to the inclusion 
of a chiral quark field in the gauge theory. Anomaly 
cancellation requires quarks to occur in vector like pairs 
so there must naturally be a partner D8. The quasi- 
normal mode spectrum we compute below will therefore 
be parity doubled in the gauge theory. 

We will again use the coordinates (10) to make infalling 
solutions regular at the horizon. 



Scalar mesons 

We first analyze the scalar mode corresponding to a ge- 
ometric fluctuation of the D8 embedding. The DBI La- 
grangian for this is 

We will parameterize the fluctuation as r = To + 
cj)(x)e~' lujh+ik ' Xa representing a mesonic excitation with 
spatial momentum k relative to the plasma rest frame 
with zero S4 spin. The Lagrangian is (a dot indicates an 
/i-derivative and a prime an x-derivative) 

C = x Wl + x 3 f(x)<p' 2 + 2#' - <£ 2 (14) 



Expanding the square root to quadratic order the equa- 
tion of motion is 




Using the regular Taylor series as initial condition we 
shoot out from the horizon and requiring our solution to 
vanish as x — * 00 we can find the quasinormal frequen- 
cies. They are shown in Fig. 2. 



Vector mesons - transverse 

We obtain the quasinormal spectrum for a Maxwell field 
on the D8 brane worldvolume, which is dual to the quasi- 
normal spectrum of vector mesons. We use the ansatz 
A^ = £,iiA(x)e lk ' X4 which is a Lorentz vector with zero 
54 spin. The equation of motion is 

da (e-*V=gF ab ) = (16) 

Fixing the gauge k^^n = one obtains the equation of 
motion for a transversely-polarized Lorentz vector. With 
this ansatz the only nontrivial equation is (for example 
choosing fc2 nonzero and the vector in the 1-direction) for 
b = 1, ie <9 Q (e-^y=5F Ql ) = 0, giving 

(xifA')'-ioj Ux-^A' - ^x-'iA s j + {uj 2 - k 2 ) x~*A = 

(17) 

Using the regular Taylor series as initial condition we 
shoot out from the horizon and requiring our solution to 
vanish as x — > 00 we can find the quasinormal frequen- 
cies. They are displayed in Fig. 3. 



Vector mesons - longitudinal 

In the zero-temperature case the transverse mesons ex- 
haust the vector spectrum. For finite temperature how- 
ever one can also identify a purely electric longitudinal 
solution of the supergravity dual Maxwell field. The lon- 
gitudinal electric field we deal with is E\ = kAo + loA\ 
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-4 1 • 

-5 



n 




1 


± 1.6016 - 0.8957 i 


2 


± 2.9217 - 1.6512 i 
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± 4.2265 - 2.4054 i 


4 


± 5.5269 - 3.1611 i 


5 


± 6.8264 - 3.9217 i 
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FIG. 3: The lowest five static vector meson quasinormal frequencies in the complex u plane. 



{k points along the 1-axis making the vector potential 
curl- free). 

There are three relevant equations of motion: first 
d a ( e ^\/~ gF ax ) = 0. Writing this out for our choice 
of metric one obtains 



and 



iojA' + ikg xl g xx A\ 



-9 U 9 x0 (k 2 A 



LukAi) (18) 



The second equation is d a (e ^^/ zr g~F a0 ) — 0. Writing 
this out for our choice of metric one obtains 




iker* J=g~g^ g Qx A\ 



Finally from d a (e 

= 



gg 1 ^ 00 (k 2 A + kwA x ) 



one obtains 



(19) 



gg^g 11 ^) - itoe^^g^g 11 ^ 
{e^^gg^g 11 {iuAt + ikA ))' 
ujkAo) 



-gg™g^ (^A, 



(20) 

The trick is to form a second order ODE for the gauge 
invariant combination E\ = kAo + wAy. We do this by 
putting the differential equations into such form as the 
coefficients of Aq and A'{ in the second and third equa- 
tions are unity then adding k times the second equation 
to uj times the third equation. We patch up the first 
derivative terms by adding zero in the form given by (18). 

(iujA' + ikg lx g xx M x + g n g x0 (k 2 A + ukA x )) = 

(21) 

We add this term with coefficient such that the first 
derivative terms add up to a multiple of kAo + ujA\ . 



The equation we finally obtain is 



Here one has 

h 



j2_ 

2x 



E'l + f 1 E[ + f 2 E x = 



,k 2 



(22) 



(23) 



£ k_ I iuj I LJ 

J2 — :c3 "T 2a; 4 / x^f 



l(uJ 2 -k 2 f) 



(24) 



After all this work, it is easy to check that this equation 
has the same k — > limit as the transverse mode - the 
quasinormal frequencies are degenerate in this limit. We 
also note that by rescaling u> —>■ X 2 uj and k — >■ Xk and 
taking A — * one finds that a normalizable, regular so- 
lution to the longitudinal equation exists for k = u> = 0, 
which is just E\ = x~ 2 . This additional mode is not 
present in any of the other spectra, and is related to the 
hydrodynamic behaviour of the field theory. 



FINITE SPATIAL MOMENTUM 



We can obtain an effective dispersion relation for our 
modes ie a function uj{k). This is done by solving the 
wave equations obtained in the previous sections for gen- 
eral complex k. Similar computations in the D3/D7 sys- 
tem can be found in [26j. 

Real momentum k corresponds to a state which is a trav- 
elling plane wave on the Minkowski spacetime of the dual 
field theory. Switching on a finite real k in the equation 
and using shooting we have found the behaviour of the 
quasinormal frequencies in the complex u> plane, ft is 
found that the states become more massive and more 
stable as k is increased. 

The results for the first three quasinormal modes for the 
scalar qq bound state are plotted in Fig. 4. In Fig. 5 we 
show the first three vector quasinormal modes for the 
transverse and longitudinal modes. These are degenerate 
for k = but behave differently as k is increased - the 
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-3 - 
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Re(w) 

FIG. 4: Lowest quasinormal frequencies at finite real spatial momentum k in the complex uj plane for the scalar qq bound state. 
The momentum k ranges from to 5.0 in steps of 0.5 (in the same units as ui) as one moves to the right in the plot. 
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FIG. 5: On the left we show the lowest vector qq bound state quasinormal frequencies at finite real spatial momentum k in 
the complex uo plane. The momentum ranges from to 2.0 in steps of 0.2 as one moves to the right (the lowest mode in the 
longitudinal spectrum, corresponding to the diffusion pole in the hydrodynamic limit, is excluded). For each mode the flatter 
trajectory corresponds to the transverse species and the steeper to the longitudinal. 

On the right we show the momentum dependence of the lowest quasinormal mode for the longitudinal vector meson. 
The equation of the line is y — log | + 2x showing the validity of the relation lu — —iDk 2 for small k, with D = |. 



main difference is that the longitudinal states become 
more stable but less massive relative to the transverse 
states for the same k as k is increased. 

Finally we note that, as mentioned above, there is an 
additional quasinormal mode for the longitudinal elec- 
tric field component of the Maxwell field on the flavour 
brane. For small k this lies close to the origin and on the 



imaginary axis. As shown in for example 25J the diffu- 
sion coefficient D for flavoured fundamental matter can 
be computed from this state. For small spatial momen- 
tum k it obeys the 'hydrodynamic' relation uj — —iDk 2 . 
Here we test whether this relation can be obtained using 
our ingoing coordinates. The calculation we do is then 
the generalization of the calculation done for the vec- 
tor meson, including a nonzero spatial momentum, and 
examining the gauge-invariant longitudinal electric field 
component. In Fig. 5 we plot the position of this pole on 



the imaginary axis as a function of k. We extract D = |. 

We note our result for the diffusion coefficient is related 
to the value obtained in (25j which is found to be D = 

77-757 . We are measuring in units of\/ — so we obtain the 

2-iTi ° y ut 

numerical value 



2tt 3 3 



(25) 



Our result therefore matches that of [25| providing a 
check on our numerics (albeit for small w and k). 



CONCLUSION 



We have found the quasinormal frequencies for a variety 
of different species (the Klein-Gordon scalar, scalar quark 



bound states and vector mesons) in the Sakai-Sugimoto 
model at high temperature. A crucial part of the analysis 
was to change coordinates so that the infalling quasi- 
normal modes become regular at the horizon so numerical 
shooting becomes straightforward. It is noteworthy that 
in these coordinates the equations for the Klein-Gordon 
scalar, the scalar qq and the transverse vector qq all have 
the form 



(x n fcj}')' - iuj (2x n ~ 3 cj)' + (n - 3)x n ~ 
+ (lu 2 - k 2 ) x n -3<f> = 



(26) 



For the Klein- Gordon scalar n = for the sigma n = 
and for the vector n = |. This means the quasinormal 
spectra look extremely similar. The only thing making 
the frequencies different is the value of n. The effect of in- 
creasing the value of n is to move the quasinormal modes 
out from the origin in the complex frequency plane. 

We also computed these states at finite real momenta 
where the modes become more massive and more stable. 
We have obtained the numerical value for the diffusion 
coefficient for fundamental flavoured matter and our re- 
sult is consistent with the calculation of 25]. 
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